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SUMMART 


The flow around cones without axial symmetry and moving at super- 
sonic velocity is analyzed. Singular points are shown to exist in the 
flow around the cone if no axial symmetry exists. The results of the 
analysis are applied to the determination of flow around circular cones 
at small angles of attack. The concept of a vortical layer around the 
cone at small angles of attack is introduced* and the correct values of 
the first-order terms of "the velocity components are determined. 

The method used is applied to cones at finite angles of attack* and 
it is shown that good agreement with experimental results can be- obtained 
from the first-order theory if the complete equation for the pressure 
distribution is used. 


INTRODUCTION 


The flow around a cone having a circular cross section and moving at 
supersonic speed has been determined by means of the assumption of small 
disturbances or by means of more rigorous methods that consider the 
existence of the shocks. The. latter methods can be applied for any Mach 
number larger than unity and have been developed by several authors* at 
first by assuming all the flow as potential flow (references 1 and 2) and 
later by also considering the variation of entropy due to the change in 
angle of attack (reference 3). By means of the development given in 
reference 3 , values of flow properties around circular cones at an angle 
of attack have been tabulated in reference 1;. The method has been 
extended in reference 3 to larger angles of attack. 

In the method given in references 3 , U* and 3 , the flow properties 
were considered continuous and were developed in Fourier series in terms 
of the angle of attackj however* the existence of a singular point at the 
surface of the cone was neglected. The derivatives of the flow properties 
were obtained by differentiating the Fourier series term by term* and the 
tezms of the series that represent the derivatives were assumed to be of 
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the same order as the corresponding terms of the integral quantities. 

For this reason, an erroneous distribution of the entropy at the surface 
of the cone was obtained. 

In this paper, the flow around the cone in the general case is dis- 
cussed, the existence of singular points in the flow is proved, and a 
different procedure for determining the flow around cones at small, but 
finite, angles of attack is developed. This procedure shows the way in 
which the values tabulated in reference U can be used if a simple cor- 
rection is introduced. The values obtained in this way are compared 
with experimental results at several values of angle of attack. 


SYMBOLS 


r,i|r,0 polar coordinates (see fig. 1) 

v r polar velocity component in radial direction (along r), 

referred to limiting velocity (see fig. l) 

v n polar velocity component normal to v r in meridian 

plane 0 = Constant, referred to limiting velocity 
(fig- 1) 

w polar velocity component normal to meridian plane 

0 = Constant, referred to limiting velocity (fig. 1) 

t time 


P 

P 

S 

7 

R — Cp 

°P 

c v 


pressure 

density- 

entropy 

ratio of specific heats (cp/cy) 


- Cy 


specific heat at constant pressure 
specific heat at constant volume 
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a speed of sound 

L projection of streamline on sphere r = Constant with 

center at center of the polar coordinate system 

V local velocity 

undisturbed velocity, referred to limiting velocity 
Vj limiting velocity (velocity for expansion in the vacuum) 

a semiapex angle of conical shock 

6 inclination of axis of conical shock with respect to free- 

stream direction 

T) inclination of axis of conical shock with respect to axis 

of body 

a angle of attack 

Subscripts : 

1 stream conditions 

a zero-order terms of Fourier series (part independent from 

angle of attack) 

b first-order terms of Fourier series (part proportional to 

angle of attack) 

c higher-order terms of Fourier series or quantities at 

surface of cone 

e quantities at external surface of vortical layer 

s quantities for polar coordinate system having axis 

coincident with axis of conical shock 

A prime is used to designate the terms of zero order in the power 
series in A0 for the quantities in the neighborhood of the meridian 
plane 0 = ji; two primes are used to designate the factor of the term 
containing A02/2 in the same power series. 



k 


NACA TN 2236 


THE FLCW FIELD FOR CONICAL FLOW WITHOUT AXIAL SYMMETRY 


In order to analyze the flow field for conical flow without axial 
symmetry at supersonic speeds, assume a polar coordinate system (r,\|r, 0). 
Call • v r the velocity component in the radial direction, v n the 
velocity component in the direction normal to r in the meridian 
plane 0 = Constant, and w the component normal to the meridian plane 
(fig. l); that is. 


v r = 


_ dr 
dt 


If the flow is conical. 


r d\Jc 
dt 


r' d8 . , 

W = — rr— Sin 

dt T 


6v r 

0 

11 

0 

6w _ 

6r 

6r 

ap = 

0 

^=0 

0S _ 

6r 


6r 

6r 


For these conditions Euler's equations become 


dv r 


w 


SlL - T 2 . y? . 0 


•n 


d\|r sin -vjr 60 


(la) 


T + + 1 3E + v v w 2 cot t = 0 

n 6t sin t 90 P 3^ 


(lb) 


w 6w 


'n 


6w 


6\jr sin >jr 60 p sin ^ 60 


+ 


v r w + vnw cot \jc = 0 


(lc) 
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and the continuity equation becomes 


0-y 

2pv r sin t + v n sin >Jr — + p sin \|r — + v n pcos\!r + w — + p — = 0 (2) 

3t 3t '30 30 


Because the energy in the flow is constant, the following, relations must 
apply: 


wv r 


y (i 

3p _ 

P 3p' 

| . L 

3v r 


3v n 


\ 

3w 

y - i\p 

30 

p 2 30y 

1 r 

30 

v n 

30 

w 

30; 

7 /I 

_ 

JP_ 3p' 

l = -L 

3vr + 


a^n . 


t 

3w 

y - i\p 

3ty 

p2 Sty ; 

1 V 

3>jr 

v n 

Sty 

w 

at; 

; equations 

(1), 

(2), and (3) results 

in 



T n 2 + w 2^ 

+ TT_ 

f* A+ i|/ 

+ ^L 

v n 2 ^ 

I **” 

3w 


L 

a2 ) 

v n 

OUU U/ 

V 

^ \ 



sin ty 

3©V 

f av n ( 

3w\ 

- n 



. 




,sin ijr 30 

aty/ = 

- u 








(3a) 


(3b) 


0 *) 


The entropy at any point of the flow can be expressed as 


s = 




+ Constant 


where p and p are the local quantities and p^ and 'p^ the stream 
quantities. Therefore, 


y - 1 3S = 13pl3p 
R 30 p 30 “ 7 p 30 


y - 1 ss = 1 ap i 3 p 

R Sty p 3\|r p 3ty 


( 5 ) 
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Combining equations (3>) with equations (1), (2), and (3) results in the 
following expressions: 






If - T n si * + Ip - ir ‘ Tn w + TrW sin * + Tn ” 008 + 


3S ^ v r 3w . w ^ v : 


yR 3\jr 


= -V T 


3>jr 


- w 


9\jc sin 30 


— + v r v n - w^cot \jr 


(6) 


Equations (6) c ombin ed with equation (la) give 


, . as as 

v n sin \lr — = -w — 
30 


(7) 


Equation (7) is general for any conical flow and defines the lines of 
constant entropy, which correspond to the streamlines. In fact, if L 
is the streamline projection on the sphere given by r = Constant, 


dS = 3S d£ + 3S d0 
dL 3^ dL 30 dL 


and, from equation (7), 


d0\ = w 

,at/L T n sin i|r 


( 8 ) 


At the surface of any conical body the component of velocity in the 
direction normal to the surface is zero and the stream moves tangentially 
to the bodyj therefore, the entropy at the surface of the body must be 
constant or must change in a discontinuous way (in which case equations (7) 
and (8) are not valid). 


THE PROPERTIES OF CONICAL FLOW WITHOUT AXIAL SYMMETRY 


In order to analyze the properties of conical flow without axial 
symmetry, consider first a polar coordinate system having its axis coin- 
cident with the axis of a circular cone at a small angle of attack 
(fig. 2) and assume that the direction of the undisturbed velocity V]_ is 
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in the plane 0 = 0, 0 = n. In this- case, the plane 0=0, 0 = n is 

a plane of symmetry of the flow and, in this plane. 


w = 0 

T n^ 0 


and 



3v, 


n 


30 


as 

30 


0 


Therefore, equation (7) shows that in the plane 0 = 0, 0 = n the • 

entropy is constant. At the surface of the cone (\|r = ijr c ), the normal 
component v n is zero and w / Oj therefore, equation (7) shows that 
the entropy remains constant also along the surface of the cone (\|r = \|r c ). 
Only at points A and B, (defined by 0 = 0, 0 = n, and \jr = \|r c ) v n -» 0 

and vi ->0} therefore, equation (7) is indeterminate. Because the body 
is at an angle of attack, the axis of .conical shock does not remain 
coincident with the direction of the Telocity V]_, and the entropy in 
the plane 0 = 0 must be different from the entropy in the plane 0 = n. 
The entropy at the cone surface therefore must be different from the 
entropy at the plane 0 = 0, from the entropy at the plane 0 = n, or 
from both. In this case-, where w = 0 and T n = 0, a discontinuity of 
entropy must exist either at A or B or both points. 


In order to find a relation between the value of the entropy at the 
surface of the cone and the values of the entropy in the meridian 
plane 0 = 0, 0 = n, the following considerations can be used: In the 

meridian plane (0 = 0, 0 = n), w, and — ^ are 2 ero because the 

30 30 

plane is a plane of symmetry of the flow field and, from equation (7), 

gg 

— = 0. Therefore, in the plane of symmetry in the zone outside of the 
singular points A and B, 


Q 

3 /3 2 s\ = 30 a 2 S 

MaO 2 / ~ v n sin t 392 


(9) 


For the case considered, the velocity component vn is negative at 
the shock or at the Mach cone and remains negative throughout the field, 
until- it becomes zero at the surface of the cone; therefore, the value 



tends to increase in absolute value as \|r 


decreases from the 


of 
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value corresponding to the value at the surface of the shock to the 

value at the surface of the cone when — is negative and tends to 

30 

decrease when — is positive. 

30 


Because the entropy r emain s constant along each streamline, the 

a 2 s 

decrease of the absolute value of — — as \jr decreases corresponds to 

a 02 

a departure of projection of the streamlines on the sphere r = Constant 


a /a 2 s] 


a 2 s 


from, the plane of symmetry, but, if j is of sign opposite to — 

a 0 Q 


ar 


the projection of the streamlines on the sphere r = Constant tends to 
converge toward the plane of symmetry as \{r decreases from the value at 
the shock to the value at the surface of the cone. 


Now, with the convention used, in figure 2, the component w is 
negative throughout the field and is zero at 0 = 0 and 0 = n. 


Therefore, 


3w 

a© 


is negative in the zone 


0 = 0 but is positive in the 


zone 0 = n, and the streamline projection tends to converge toward the 
zone of point A and diverge from the zone of point B. Because of the 
departure of the streamlines from the plane 0 = n, the entropy in 3 one B 
remains constant, and, therefore, the entropy at the surface of the cone 
is equal to the entropy at the meridian plane 0 = nj at point A a dis- 
continuity of entropy exists from the value corresponding to the plane 
0 = 0 to the larger value existing in the plane 0 = n.^ All the pro- 
jections of the streamlines converge at point A where th*e entropy is not 
single-valued. Because v n approaches zero near the cone, equation (8) 
shows that all the streamline projections tend to become parallel to the 
line \(r - Constant in the zone near the cone and converge at A. The 
value of v n in the meridian plane 0 = 0 can change sign and can be 
positive in the neighborhood of the point A (case of large angles of 
attack). In this case the right-hand side of equation (9) changes sign 
and the singular point moves away from A in the meridian plane 0 = 0 
and occurs at the other point of the meridian plane where v n is also zero. 
(At the shock or at the Mach wave v n is negative; therefore, another 
singular point where v n = 0 must exist. ) 


It is interesting to observe that singular points must exist in any 
supersonic conical flow without axial s y m m etry. Considerations similar 
to those used for the case of circular cones at arf angle of attack can 
be extended to other cases, and it can be shown that the streamlines that 
are tangent to the body start from points of ‘the shock and meet the body 


/ 
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at points where the component of velocity perpendicular to the radius 
and tangent to the body vanishes and has a positive derivative in the 

streamline direction ^equivalent to the condition of positive Con- 


vengency of streamlines occurs and, therefore, the points are sing ular at 
the points where this component vanishes and has a negative derivative in 
the direction tangent to the body, while the component normal to the body 
also vanishes and has a negative derivative in the direction normal to 
/ 3v n \ 

the body (equivalent to v n = 0 and — — negative], 

\ 3 | / 


For example, the conical body of figure 3 has two planes of symmetry, 

3w 


AA' and BB' when w and v, 


n are zero, but at the points BB*, — is 

30 


positive, while at AA* , it is negative. Therefore, because — — is 

3\|r 

negative at AA' and BB 1 , the points AA* are singular points and the 
entropy at the surface of the body is determined by the shock strength at 
the points CC ' . 


DETERMINATION OF THE FLOW AROUND CIRCULAR CONES 
AT SMALL ANGLES OF ATTACK 


In order to determine the flow around circular cones at small angles 
of attack, consider a polar coordinate system, the axis of which is coin- 
cident with the body axis. At the surface of the body the velocity com- 
ponent v n is zero and in the neighborhood of the body is very small j 
therefore, the terms v n 2 /a 2 can be neglected with respect to unity. 

If the angle of attack is small, the component w is also small 
and the terms w^/a 2 can also be neglected. On the basis of this 
approximation, in the neighborhood of the body, equation (lj.) can be 
expressed as 


—5 +■ — a g ... . . = 0 

3\jr sin >Jr 30 


2v r + v n cot i|r + 


( 10 ) 
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This equation permits a particular solution of 
physical considerations 

v r = v ra + v^ cos 0 +]T v r c 

1 

Vn = cos 0 +£jm c 

w = wjj sin 0 + Y~ w c sin m0 


this type chosen from 
cos m0 1 


cos m0 > 


( 11 ) 


where v r , T n , and W5 c are functions only of \|r, are 

constant for constant values of ty, and must be chosen in a form that 
satisfies the boundaiy conditions. 


Consider now a conical shock havin g circular cross section and semi- 
apex an g le ct (fig. U) • Consider the cone inclined at an angle 6 with 
respect to the undisturbed velocity with a polar coordinate system, the 
axis of which is coincident with the axis of the conical shock. If v rs . 


v„ , and w„ 
n s s 


are the velocity components referred to the limiting 


velocity V l in the new coordinate system (r s , 0 S ), from the shock- 
wave relations the following equations result: 


v rs = cos 6 cos cr + Vi sin 6 sin ct cos 0 S 


w s = -Vi sin 6 sin 0 S 


— Vyi— = 


7 - 1 . 


1 " v r s 2 - w s 2 


n s y + 1 Vi cos 6 sin ct - Vi sin 6 cos ct cos 0 S 

% 

If 6 is assumed to be small and terms of the order of 6^ are 
neglected, these equations become 


Vr s = Vi cos ct + 6V1 sin ct cos 0 s 
w s = -V16 sin 0 S 


7-1/^ 

L - Vi^cos^d 

w - 1 

7 + 1\ 

Vi sin ct 

/ 7 + 1 


/ 1 - Vi^cos^ctN 

6 cos ct|-2Vi- + - /cos 0 S 

y Vi sin 2 CT / 


> ( 12 ) 
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The flow behind a circular conical shock inclined at a small angle 6 
with respect to the undisturbed stream can therefore be expressed in the 
form 


T r s - Tr as + Bv^ 00 s 0 S 

T “b s ' T % s + 6T “b a 008 08 
w = 6wb g sin 0 S 


> (13) 


where all the terms containing v^, v^, and w^ are small, that is, 

of the order of 6. 


If the axis of the conical coordinates is rotated at an angle t) 
of the same order as the angle 6 and terms of the order of 6* are 
neglected, the velocity components v r , v n , and w referred to the new 
axis (fig. U) become 


3v- 


r a e 


vr = v r + ■ ° T) + Sv^ cos 0 

a s 3\Jr s D s , 


dV; 


Vn ' Tna s + ~5fT T ’ + 6Tnb8 003 0 


(U.) 


w => 6wb g sin 0 

where 0, the coordinate referred to the new axis,* is given by 

cos 0 = cos 0 S - T) cot i|r s sin^0 s 


Equations (llj.) show that the flow behind a circular conical shock 
inclined at a small, angle 6 with respect to the direction of the undis- 
turbed stream and at an angle r] with respect to a conical body can also 
be expressed in the form given in equations (11), where, if the angles T] 
and 6 are s mall and the terms of the order of tj 2 and 6^ are neg- 
lected, only the terms having the subscripts a and b must be con- 
sidered. Therefore, equations (11), when the terms with the subscript c 
are neglected, are valid for small angles of attack (a * 6 - T]), and a 
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conical circular shock is consistent with the solution chosen for the 
flow around the 'body. At the surface of the cone the assumption 
that w^ is small corresponds to the assumption that only the first- 
order terms of angle of attack are considered. The conical shock is 
inclined at an ang le T} with respect to the circular cone. 

This analysis is similar to the analysis of references 2 and 3., No 
assumption, however, has been introduced for the entropy distribution^ 
only the velocity components have been considered to be in the form of 
equations (11), and no limiting assumption has been introduced for the 
derivatives. In reference 3, in addition to the equations 


the expressions 


Vf, = V-r + OV- COS 0 
x a 

(l5a) 

v n = X + OT nb 003 9 

(15b) 

w = aw^ sin 0 

(15c) 

P s P a + ^Pb cos 9 

(l6a) 

P 3 P a + cos 0 

(l6b) 

S = S a + aS^ cos 0 

( 16 c) 


have been used, and the derivatives of entropy, pressure, and density 
have been obtained from differentiation of equations (16). In this way 
a solution has been found which gives values of entropy variables along 
the cone surface and constant in each meridian plane, while the entropy 
actually remains constant along the cone surface and changes in the 
meridian plane. An incorrect entropy distribution has therefore been 
obtained at the surface of the body. 

In order to analyze in more detail the significance of equations (l6) 
and their inconsistency with the approximation "considered in references 2 
and 3, consider equation (7). In the plane of symmetry 0 = 0 or 0 = n, 

— = Oj therefore, S a + aS-^ 'or S a - aS^ of equation (16c) must' remain 
3\jr 

constant. 
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form 


Consider now the plane 0 = n and express the entropy S in the 


S = S' - S" 


A0 2 


(17) 


which satisfies the condition of symmetry. Because of equation (7), 

dS* 


d\Jr 


and, from equations ( 16 ) and (17), 


= 0 


(18) 


S» = S a - otSb (19) 

and 

oS b = -S» (20) 


From equations (7) and (17), 

v„ 3in JML . SSL i£j. rs« Ae 

\d>jr d \| r 2 ) 

as> 

However, — — = 0: therefore, 

A0 SS^ _ wS" 

2 ity v n sin y 

or, since from equation (15c) w = -awb A0, 

8S a 2awbS» 

9^ v n sin \|r 

By use of equation (20), in the neighborhood of the plane 0 = n, 

l 

a = 2w b s b° 2 = 2wS" 

v n sin t|t v n sin ^ 


0 
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In reference 3, the term 


wS" 


v n sin t 


to be of the order of 
8S 


has been considered everywhere 

9% 


sr and has been neglected; hence, = 0 

d\|r 


and 


at 


= 0. However, neglecting this term is correct only when the 


ratio wfc/v n is of an order different from 1/a and, therefore, 
when M » 0. Near the surface of the cone, v n approaches zero, 

and, therefore, the term becomes large and cannot be neglected. 

at 

T as n 

The extent of the field around the cone where the term is of the 

at 

order of a can be easily determined. 

Consider a polar coordinate system having its axis coincident with 
the axis of the cone. At the surface of the body v n is zero; therefore, 
in the neighborhood of the body the velocity component v n can be 
expressed in the form 


7 n = (—) ^ 

\9t /tc 


or, ty use of equation (10), in the form 






2'v r 


3w 
sin t 96. 


•) Ut)v 


( 21 ) 


Therefore, v n is of the order of a when (At)y is the order of a. 

* as 11 

In this conical layer of thickness (At)v of the order of a, — — is 

3t 

also of the order of a or larger. In this layer, which can be called 
the vortical layer, the term 

from equations (18) and (19) 


3S" 



In order to investigate th*e effect of this vortical layer on the 
velocity and pressure distribution at the surface of the cone, consider 


as 


at 

as 


— also is of the order of a because 


at 


— can be" shown to be of the same order 
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equations (2), (1;), (6), and (7). If the velocity components in the 
neighborhood of the meridian plane (0 = n) are expressed as 



the following expressions can be obtained: 

as» 2w"S" ' 

■ ■ ■ ■ 3 — ■ 

3t v n * sin t 



(22a) 

(22b) 

(22c) 

(23) 

(210 


aw» 

at 


" V sin + Vn ' Vn " + W ” Vr ' S±C1 + Vn ' W " CO ? ♦ + “£f“) (2 ^ 


Qvn” 1 

at 


1 - 


v n » 


l«2 


,vn Ah ,v n ' 1 , v n ,2 \ „ , . v n'^ 

l~ - A ^r) - v "" cot t - - -jrj - 


w 11 L + 2w n ^ _ ^n'^n” + 2v n 1 sin t 3wA 
sin t\ a*2 a'2 a'2 3t / 


,i2 


. (26) 


where 


A = 


-7 l-nr II — V T'or II + wll2 

r r n n 

ITi 



2 


(27) 
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v r » + v n ' cot -vjr + 


w" 


B = - 


sin \| r 


1 - 


/ 


a 


,2 


and 



a «2 1 5S a 

vr a y R ^ 


(28)' 


(29) 


Therefore, all the derivatives of the velocity components of zero 
and first order are affected by the entropy variation in the meridian 


. as 


3S" 


plane. If the terms — - and are of the order of 


at 


at 


a, the cor- 


O 

rection required for the entropy variation is of the order of aA because 
the thickness of the vortical layer is of the order of a and, 'therefore, 
must be considered only if higher-order terms are considered, as in refer- 

dS 

ence S>. A correction is also required in this approximation when ^ 

t&C5|| 

and — — become of zero order. 

at 


The only place where the terms 


are of zero order is 


OSl. and SSI 

3* a* „ 

is near the surface of the cone in a layer of thickness aA where v n 
tends to the order of a 2 (equation (9)). Because this layer is of the 
order of a 2 , it can be considered in this approximation, which neglects - 
terns of the order of a 2 or higher, that across this layer the pressure 
distribution remains constant, but an abrupt variation of entropy occurs j 
therefore, in this approximation the phenomenon can be represented as in 
references 3 and U, where the entropy remains constant in every meridian 
plane -until a vortical layer of infinitesimal thickness is reached at the 
surface of the cone across which a variation of entropy occurs from the 
value S a + aSb cos 0 to the value S a - aSb that exists at the surface 
of the cone. Across the layer a variation of density and velocity com- 
ponents occurs and can be easily determined. 


Gall v r , w e , S e , p e , p e , and a e the quantities at the 
© 

external surface of the layer (these are the quantities tabulated in 

reference I 4 .) , and v r , w c , S c , p c , and p c the quantities at the 

c 

surface of the cone. Because it has been shown that 


Pe “ Pc 


(30) 
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then 


Vr 2 ♦ we 2 - v r 2 - i_ a „2 h—fs. 


y - 1 ' cp 

where Se - Sc is the entropy jump across the layer. 
Now, 

v r = v ra + OVr^ cos 0 


(31) 


and in the plane 0 ** n 


w = awb sin 0 


T r e = v r c 


w e = w c = 0 


Se ~ S c 


or 


( T r a - onr i’ b ) e - ( T r a - ^rb) c = T r’ 

(s a - aSb) e = (s a - aSb) c 53 S' 


(32) 


Therefore, if terms of the order of are neglected, equation (31) 
becomes 


Tr * ( Tr b)e * T r'( T rb) 


a 2 fb 

C 7-1 Cp 


(33) 


or 


( Tp b)c - ( Vr b)e 


a 2 S b 
(7 - l)v r * c p 


(3U) 
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However, from, equation (la), at the surface of the 

3v r 

w = — — — — 
sin \jr 30 

or 


b sin ijc c 

therefore, the values of v r and can he determined from equa- 
tion (3U) and from the tables of reference 1* where 

Sb ' ° v (iv ' 7 (36) 

and 

a e 2 = " Vr ‘ 2 ) (37) 

The values v ra , v-^, and having been determined, the values of v r 

and w can be obtained at any part of the cone from equations (l5) . 

* 

From the velocity components, the correct pressure distribution can 
be determined since the entropy remains constant and is equal to S a - aSfo 
at the surface of the cone. 


cone. 


(35) 


THE NUMERICAL DETERMINATION OF THE FLOW FIELD AROUND 
CONES AT SMALL ANGLES OF ATTACK 


The method presented permits the determination of the pressure dis- 
tribution around the cones with the assumption that the terms of the 
order of a? can be neglected. The pressure at any point can be obtained 
from the equation 


_P_ 

Pi 



1 AS 


( 36 ) 



NACA TN 2236 


19 


where ¥ is the local velocity corresponding to the pressure p 
and AS is the increase of entropy with respect to the stream condi- 
tions where pi and V]_ exist. From equations (15>)> V 2 is found 
to be 


V 2 = 


+ a 2 v r ^ 2 cos 2 0 + 2crv ra v r ^ cos 0 + a 2 Wb 2 sin 0 


(39) 


However, in equation (39) the terms of the order of a 2 having the form 
2v r a v r c cos . m0 have been neglected (see equation (11)). If all the terms 

of the order of a 2 are neglected, V 2 becomes 

Y 2 = v r a 2 + 2^^^ cos 0 (1*0) 


Equations (39) and (1;0) are different in terms of the order of a 2 ; 
however, for finite angles of attack good agreement is obtained only if* 
equation (39) in which some of the a 2 terms are retained is used in 
equation (38). 


The reason for the better approximation given by equation (39) 
be understood if the magnitude of the terms containing a2 in the 
expression of V 2 i 3 considered. Along the surface of the cone, awfc 
is given by 


awb = 


OT rb 

sim|r c 


(ia) 


For a finite value of a and a small value of >|r c , a2w^2 i s of the 

same order as av^ because sin 2 \|r c is also small; therefore, a 2 wb 2 

can have an effect on the velocity and pressure distribution of the same 
order as the term av Ts r r ^ which is the only term retained in equa- 
tion (U0). . Because the term a^w 2 is correct and significant, it can 
still also be retained if other terms in a 2 are neglected, and, there- 
fore, equation (39) is the expression that must be used for finite angle 

of attack. (For example, for a 10° cone (i|r c = 10°) sin 2 >|r c = — 

CL ^ \ 

and „ = 1 for a = 1.7S>°, which can be considered a small angle.) 

sinr^ J 


In reference 3, equation (16a) has been used in the derivation of 
the method; however, the use of this equation is not necessary, as is 
briefly shown in the following paragraph. 
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Consider a conical circular .shock having the semiapex angle a and 
incline d at an angle 6 with respect to the undisturbed stream. In the 
neighborhood of the plane 9 - n the velocity components can be expressed 
as in equations (22) and at the shock are given by (from equations (12 JJ 


v r t = Vi cos a - V]6 sin a 


v r " = _6Vi sin a 


w" = -V 16 


m 


(1*3) 

m 


- V -, 2 


1 - V 


T n' = 


y - 1 

7 + 1 


COS^CT y - 1 

+ - 6 cos a 


sin o 


y + 1 


( 1 - ¥-i 2 cos 2 cn 

(- 2 ?! + P )( h$) 

\ x ?! sin 2 a / 


1 / 1 - ¥i 2 cos 2 o\ 

v » = L 6 cos al-2Vn + — — 7-5 ) 

n y + 1 \ ¥1 sin 2 cr / 

and S' and S" are given by 

S' -.oy log e jl + Mi 2 + l) [y 1 2 sin 2 (a + 6) - v n ' 2 ] 

— Vji* 

(°p ' c ^) loge Vl sin(o + 6) 


0 * 6 ) 




0*7) 


S" = 26- 


¥-i 2 cos a sin a - v n ' —• Ah"/ 6 cos ct\ 

£ * (cp ' v" ’ 

— 1 TUT- 2 ' 


0*8) 


1 + Z-=_A Ml 2 


an the velocity components are referred to the limiting velocity, and 
Y 1 i s the undisturbed velocity also referred to the limiting velocity. 

In the meridian plane 0 = n ' the entropy S' is constant; therefore 
when w = 0, from equation (la.). 


3v r » 

cty 


s TT I 

v n 


(1*9) 
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If equation (U) is applied to the meridian plane 0 = n and equa- 
tion (21) is used, the radius of the hodograph diagram at any value of \jr 
smaller than a is given by (reference 6) 


where 




1-— i 

a*2 / 



= 1 - 




(50) 


(5i) 


Therefore, at the point ijr - A\Jr in the meridian plane 0 = it. 


( v r' )i{f-A>{f = “(Tn'Xjr sin A\(r + (v r l - R)^ cos Af +'(R) 1 j f (£2) 

and 

(v n * )\Jr— A^Jr " (vn* )y cos Aijc + (v r ‘ - R)^ sin At|t (53) 


From equations (50) to (53), the velocity components v r * and v n * 
can be determined if the component w n is known in the- meridian 
plane 0 = n. Since — — - is given by equation (25), the value of w" 

OY 

can be determined for any value of ijr. But v Pa and v^ are the 
quantities obtained for zero angle of attack at the same coordinate \|r 
of a coordinate system in axis with the conical shock, and 

v r " = v r » - v r& 

▼n" = T n l " 

therefore, the entire flow field can be determined until v n ' becomes 0. 
The value of \J t for which v n ' = 0 corresponds to % + a - 6 and, 
therefore, gives the value of a. 

Considerations similar to those used for cones can be used for the 
characteristics method presented in reference 6. In this case the 
pressure at the surface of the body can be obtained from the complete 
equations, and the vortical layer must be considered in order to obtain 
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the correct distribution of entropy. The application is the same because 
the entropy does not change along each streamline. 


COMPARISON WITH EXPERIMENTAL RESULTS 


In order to have an indication of the accuracy that can be expected 
from the first-order theory, theoretical results have been compared with 
seme experimental results available. The theoretical results have been 
obtained by using the values of reference U for the conditions outside 
the vortical layer and the pressure at the surface of the cone has been 
determined in the following way: 


From tables of reference U the value of 6/a has been determined 
(6/e of tables). The position ijr in the plane 0 = n of the conical 
body in the shock coordinate system is 


t = tc + a ~ 6 


m 


The value of v r „ at \}r c + a - 6 has been determined from the tables 


of reference k ^v r at \|r is given by us - 2us }. The value 

of v r " has been obtained from the tables in shock coordinates 


(a - 6)2j 




(55) 


Because — - =0 at \j r = \|r c , the value of J^v-q^) JJ ^ can be taken equal 
to the value of vp^ at \|r and, therefore, to -x at of the 
tables. Then ( v r a ) c > ( v rb)c J (^Jc 1x1 * iie ^ 0< ^ r coordinate system 
have been obtained by means of equations (32), (3U), and (35) • 


From equation (36) St, has been obtained (it has a negative value), 
and from tables of reference 7 S a can be determined - for example, 
from the value of the angle of the shock obtained from reference 7 - and 


AS = Sa - aSfc 


( 56 ) 


The pressure has been obtained from equation (38). 
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In figure 5 a comparison is presented for a cone of \|r c = 7.5° at 
M = 1.6 and four angles of attack. The experimental data are obtained 
from tests performed in the Langley i|- by i;-f oot supersonic tunnel. 

For comparison the values given by references 1* and 5 and by linear 
theory are also shown. In figure 6 a comparison is presented for 
a cone of i|r c “ 10° at M =-6.86 and two angles of attack. The experi- 
mental data have been obtained from tests performed at the Langley 
11-inch hypersonic tunnel. The agreement in both cases is good, even at 
angles of attack where it would be expected that higher-order ter ms would 
be important. 


CONCLUSIONS 


The flow around cones without axial symmetry at supersonic velocity 
has been analyzed. Singular points which complicate the analysis of the 
flow field were shown to exist in the flow. The results of the analysis 
were applied to the determination of the flow around circular cones at 
an angle of attack. The concept of a vortical layer around the cone at 
sm a ll angles of attack has been introduced, and the correct values of the 
first-order terms of the velocity components were determined. 

The method determined was applied to cones at finite angle of attack, 
and it is shown that good agreement with experimental results can be 
obtained from the first-order theory if the complete equation for the 
pressure distribution is used. .. The analysis can be extended to the 
application of the characteristics method around bodies of revolution 
at small angles of attack. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va., September 26, 1950 
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Figure 2 .- The singular points at the surface of the cone. 
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Figure 4.- The position of the conical shock with respect to the 

conical body. 
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(a) a = 4®. 


O Experimental data 

Theoretical data 

Linear theory 

—Theoretical data from reference 4 

Theoretical data from reference 5 


(o) «= 8 


(b) e:* 6°. 


(d) oo » 12°. 


Coordinate 8 


Coordinate 8 


Figure 5*- Comparison of experimental pressure distribution over the 
surface of a cone, with the theoretical pressure distribution. Semi- 
apex a n gle of cone, 7*5°i M = 1.6. (Experimental data obtained from 
the Langley 4- by 4-foot supersonic tunnel . ) 
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O Experimental values for two 
□ diffe^pnfc values of r 
along the oone surface 
Theoretical values 



Figure 6.- Comparison of experimental pressure distribution over the 
surface of a cone with the theoretical pressure distribution. Semi- 
apex angle of cone, 10°; M = 6.86. (Experimental values for two 
different values of r along the cone surface^ data obtained from 
the Langley 11 -inch hypersonic tunnel. ) 
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